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Abstract 

To investigate the formation and the propagation of relativistic shock waves in viscous gluon 
matter we solve the relativistic Riemann problem using a microscopic parton cascade. We 
demonstrate the transition from ideal to viscous shock waves by varying the shear viscosity to 
entropy density ratio 77/5. We show that an tj/s ratio larger than 0.2 prevents the development of 
well-defined shock waves on time scales typical for ultrarelativistic heavy-ion collisions. These 
findings are confirmed by viscous hydrodynamic calculations. 



1. Introduction 

Recently, in the context of the jet-quenching phenomena QJ], exciting jet-associated particle 
correlations [2] have been observed in heavy-ion collisions at BNL's Relativistic Heavy-Ion Col- 
lider (RHIC). This indicates the formation of shock waves in form of Mach Cones [3] induced 
by supersonic partons moving through the quark-gluon plasma (QGP) and could give a direct 
access to the equation of state of the QGR 

Shock waves can be observed if the matter, where the shock waves are induced, behaves like 
a fluid. The large measured elliptic flow coefficient V2 0] implies that the QGP created at RHIC 
could be a nearly perfect fluid with a small viscosity. Calculations of viscous hydrodynamics |5] 
and microscopic transport theory ((J have estimated the shear viscosity to the entropy density 
ratio rj/s to be less than 0.4. Still, it is not known, if this upper limit of the 77/ s ratio is sufficiently 
small to allow the formation of relativistic shock waves. 

In this work we address the question, whether and when relativistic shock waves can develop 
in viscous gluon matter for given rj/s values. For this purpose we consider the relativistic Rie- 
mann problem. Its initial conditions are two regions of equilibrated matter with different constant 
pressure separated by a membrane which is removed at t — 0. The matter evolves in form of a 
one-dimensional expansion. We solve the relativistic Riemann problem employing both the par- 
ton cascade BAMPS (Boltzmann Approach of MultiParton Scatterings) [7] and the vSHASTA 
approach (viscous SHArp and Smooth Transport Algorithm) [8] of viscous hydrodynamics. We 
demonstrate the transition from ideal shock waves to free diffusion by varying the rj/s ratio from 
zero to infinity. We estimate the upper limit of the r]/s value, for which shocks can still be 
observed experimentally on the time scale typical at RHIC. 
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2. BAMPS and vSHASTA 



BAMPS is a microscopic transport model solving the Boltzmann equation 

p u d„f(x,p)=C(x,p) (1) 

for on-shell gluons with the collision integral C(x, p). The algorithm for collisions is based on 
the stochastic interpretation of the transition rate J7[] . In this study, we consider only binary gluon 
scattering processes with an isotro pic cross section, which is adjusted locally at each time step 
to keep a constant rj/s value 

vSHASTA solves the Israel-Stewart (IS) equations of dissipative hydrodynamics. In 1 + 1 
dimensions IS equations reduce to |0] 

d,T 00 + d z (vT 0l> )=-d z (vP + vff), (2) 
d,T Qz + d z (vT Qz )=-d z (P + n), (3) 

yd,n + yvd z ii=—{n m - ft) - ^ + D In^j , (4) 

where 6 = and D = u^d^. For vanishing shear viscous pressure, n, Eqs. (0 and (0 reduce 
to the equations of ideal hydrodynamics. 



3. Development of shock waves 

Figure [T] shows the solution of the relativistic Riemann problem for various rj/s values as 
computed with BAMPS. The results demonstrate a gradual transition from ideal hydrodynamic 
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Figure 1: (Color online) The solution of the Riemann problem. At t = 0, the pressure is Pq = 5.43 GeVfm : for z < 
and P4 = 0.33 GeVfm -3 for z > 0. The left panel shows the pressure and the right panel the velocity at time t = 3.2 fm/c. 



limit (rj = 0) to free streaming of particles (77 = 00). The ideal solution II 1311 is reproduced with 
a very good precision for rj/s — 0.001. We see that a shock wave is traveling into matter with 
lower pressure (on the right) with a velocity v s h oc k larger than the velocity of sound c s . The 
region behind the shock front is the so-called shock plateau, where matter is equilibrated and 
moves collectively with v p i at , shown on the right panel of Fig. [T] Simultaneously a rarefaction 
wave is propagating with c s into matter of higher pressure (on the left). 
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A larger 77/ s value results in a finite transition layer where the quantities change smoothly 
rather than discontinuously as in the case of a perfect fluid. Furthermore a non-zero viscosity, if 
large enough, impedes the formation of a shock plateau and a clear separation of the shock front 
from the rarefaction fan. 
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Figure 2: (Color online) Same as in Fig. fj] Results are obtained using BAMPS and vSHASTA. 



In Fig.|2]we compare the results from BAMPS and vSHASTA for tj/s = 0.01 and 0.1. We 
see a perfect agreement for rj/s = 0.01, whereas for larger value of rj/s — 0.1 small deviations in 
the region of the shock front and rarefaction wave are found. The reason for the difference is that 
in these regions the local Knudsen number Kg = A^d^u 11 111411 is large and thus the applicability 
of IS equations is questionable. Transport calculations do not suffer from that drawback. 

4. Time scale of formation of shock waves 
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Figure 3: (Color online) Same as in Fig. \T\ Results are obtained using BAMPS for );/* = 0.1. 



The formation of a shock wave takes a certain amount of time, as demonstrated in Fig.[3]for 
rj/s = 0.1. At early times a shock has not yet developed, the profile looks like a free streaming 
of particles. But at later times we observe the creation of a shock plateau. Formally we define 
the time of formation of the shock plateau when the maximum of the velocity distribution v(z) 
reaches the value v p i at of the ideal-fluid solution in Fig.Q] From the right panel of Fig. [3] we see 
that this happens at t — 3.2 fm/c. 
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In Ref. [9] we have found the scaling behavior of the Riemann problem: the time scale of the 
formation of shock waves t/ is proportional to the rj/s ratio. From Fig. [3]we also infer that, for 
rjl s > 0.1, a shock plateau has not yet developed at t — 3.2 fm/c, whereas for J]/s < 0.1, it has 
already fully formed. For 77/s > 0.2 we obtain tf > 6.4 fm/c, which is the typical lifetime of the 
QGP produced in heavy-ion collisions at RHIC. Therefore, an imperfect quark gluon fluid with 
rj/s > 0.2 will prevent the observation of Mach Cone signals at RHIC. 

5. Conclusion 

Using the parton cascade BAMPS we have solved the relativistic Riemann problem. The 
transition from ideal-fluid behavior to free streaming is demonstrated. Numerical results from 
BAMPS agree well with those obtained from viscous hydrodynamical calculations using vSHASTA. 
Considering the scaling behavior we found that the formation of shock waves in gluon matter 
with rj/s > 0.2 probably takes longer than the lifetime of the QGP at RHIC. 
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